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Abstract 

The Joule heating effect on graphene electronic properties is investigated by using full-band 
Monte Carlo electron dynamics and three-dimensional heat transfer simulations self-consistently. 
A number of technologically important substrate materials are examined: Si02, SiC, hexagonal 
BN, and diamond. The results illustrate that the choice of substrate has a major impact via the 
heat conduction and surface polar phonon scattering. Particularly, it is found that the poor ther- 
mal conductivity of Si02 leads to significant Joule heating and saturation velocity degradation in 
graphene (characterized by the so-called l/y/n decay). Considering the overall characteristics, BN 
appears to compare favorably against other substrate choices for graphene in electronic applica- 
tions. 

PACS numbers: 72.80.Vp,72.10.Di,72.20.Ht,65.80.Ck 
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Superior electronic properties are one of the main attractions of graphene in device appli- 
cations. The massless Dirac Fermions originating from a linear dispersion relation imply high 
electron mobilities and drift velocities — an ideal trait for devices in high- frequency /high- 
speed operation.- Nevertheless, carrier motions are subject to scattering by perturbations 
such as lattice vibrations or impurities in most (if not, all) realistic conditions. Particularly, 
the electron interaction with the lattice is essentially a thermalization process from a sys- 
tem point of view. As electrons gain energy from an external source (such as an electrical 
bias), a part of the excess energy is transferred to the lattice via phonon emission. Subse- 
quent increase in the lattice temperature (i.e., the Joule heating) acts as a counter weight to 
limit further energy gain from the source by causing degradation in the electronic transport. 
Eventually, a balance is reached and the system approaches the steady state. Thus, the 
details of heat dissipation including the properties of its primary path (i.e., the substrate) 
could have a major influence. This is even more so in graphene based structures,- where the 
two-dimensional (2D) nature dictates a large interface with the substrate compared to the 
volume. 

In this Letter, we theoretically investigate the effect of Joule heating in graphene. Specif- 
ically, the impact of different substrates on graphene electron transport properties are ex- 
amined with four technologically important substrate materials: SiC>2, SiC, hexagonal BN 
(A-BN), and diamond. While SiC>2 is the most commonly used substrate due to its compati- 
bility with conventional technology, 3,4 SiC has seen its use in producing large-area graphene 
by solid state graphitization.™ Recently, /i-BN has drawn attention for its structural similar- 
ity to graphene — a much desired condition for high quality samples.- On the other hand, 
diamond can also be a candidate. Beside the anticipated affinity with graphene, it is one 
of the most thermally conductive materials. In the analysis, we consider a model problem, 
where a small graphene sample is placed on a 2D plane of relatively thick dielectric or sub- 
strate (300 nm), which is in turn on top of a bulk Si layer. The graphene film is subject 
to a uniform electric field, generating excess heat that must be dispersed through the layers 
underneath. To obtain the self-consistent solution across the structure, we solve simultane- 
ously the 3D heat transfer equation (including the estimated interfacial thermal resistance) 
and the Monte Carlo electron dynamics in graphene. 

The Monte Carlo simulation developed in this study takes into account the complete elec- 
tron and phonon spectra in the first Brillouin zone. Specifically, both the graphene phonon 
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dispersion and its interaction with electrons are obtained from the density functional the- 
ory calculations,- whereas a tight binding model is used for the electronic energy bands.- 
In addition, graphene electron interactions with charged impurities (5 x 10 11 cm" 2 ) on the 
substrate surface 9 - as well as the surface polar phonons (SPPs) are included.— For simplicity, 
the phonon system is assumed to reach thermal equilibrium fast enough so that the electron- 
phonon scattering rates, for both intrinsic graphene phonons and SPPs, have the tempera- 
ture dependence based on the Bose- Einstein distribution: i.e., N g = l/[exp(-£^-) — 1], where 
u p h is the phonon frequency, ks the Boltzmann constant, and T the temperature. 

The electron energy transferred to the lattice vibrational modes increases the lattice 
temperatures of graphene and the substrate. Specifically, the interaction with graphene 
phonons (e.g., emission) leads to the elevation of graphene lattice temperature (T g ), whereas 
the substrate has contributions from both the direct excitation of SPPs (i.e., electron-SPP 
scattering) and the heat conduction from the graphene lattice. As summarized above, the 
thermal part of the self-consistent model utilizes a 3D heat transfer equation in the substrate; 
V • [k(x, y, z)VT(x, y, z)] = 0, where « is the thermal conductivity of the material. The 
values used in the calculations are Ksi0 2 = 1-4 Wm~ K -1 , ksic — 370 Win K and 
K dia — 1800 Wm K for SiC>2, SiC, and diamond, respectively.— 1 ^ Unlike the first three, 
the thermal conductivity of /i-BN is anisotropic with a large difference in the in-plane and 
the out-of-plane direction due to the layered nature; Kh-BN(x,y) ~ 300 Wm _1 K _1 and 
ft-h-B^i?) ~ 2 Wm _1 K -1 .— The corresponding details on Si can be found in Ref. Q. 

In a heterogeneous system, there is an extra thermal resistance r gs at the interface between 
graphene and the substrate (i.e., the so-called Kapitza resistance). An experimental mea- 
surement reported r gs ranging from 5.6 x 10 -9 to 1.2 x 10~ 8 Km 2 W _1 in the graphene/Si02 
structure.— Interestingly, a first principles calculation conducted very recently also suggests 
similar numbers for the graphene interface with /i-BN and SiC.- 6 As such, a typical value of 
8.8 x 10~ 9 Km 2 W _1 (i.e., the median of the range observed for SiC^) is adopted in this study 
for all four substrate materials. Concerning SiC, however, the situation is more complex. 
When it is exposed to the air, hydrogen tends to be adsorbed and terminate the dangling 
bonds of Si lonely atoms in order to form a stable surface.— The surface is also passivated 
intentionally to reduce the interface states. As there is an indication that this could cause 
a drastic increase in r 9S ,~ an additional case of fully hydrogen terminated SiC (SiC-H) is 
considered with r gs = 7.9 x 10 -8 Km 2 W _1 to gauge the impact. In the real situation, the 
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surface is more likely to show partial termination (i.e., somewhere between the cases of SiC 
and SiC-H). 

Then, the total power dissipation per unit area across the interface between graphene 
and the substrate can be expressed as P t = (T g — T s )/r gs + P spp , where T s is the lattice 
temperature of the substrate at the interface and P spp is the power transferred via the direct 
SPP scattering. As Pt should be equal to the net power loss by graphene electrons in a 
steady state, both this quantity and P spp can be estimated from the Monte Carlo simulation, 
providing the necessary boundary condition for the heat transfer equation. Finally, a self- 
consistent temperature profile (including T g and T s ) is obtained by an iterative process. 

Figure 1 illustrates the potential impact of Joule heating in graphene on different substrate 
materials. Two sets of data are provided with the electron density of 1 x 10 12 cm -2 and the 
graphene sample dimension of 1 /im x 0.5 /jm. One set of results, shown in lines, represents 
the drift velocity versus electric field when Joule heating is ignored by fixing T S =T S =300 K. 
The other set, in data points, examines the same curves with the Joule heating effect taken 
into account. Of the cases under consideration, the most drastic changes appear in the 
graphene/Si02 structure. While the deviation is minor in the low-field region, the saturation 
velocity v sat degrades substantially (e.g., by about 20% to 3.9 x 10 7 cm/s at 30 kV/cm). 
In comparison, the velocity-field curves show little impact of Joule heating for all other 
substrates. Only SiC-H shows a minor influence; the rest including SiC (unpassivated; not 
shown in Fig. 1) appear virtually unaffected. Here, it is also interesting to note that the 
graphene-on-diamond structure actually has the lowest v sat independent of Joule heating. 
The origin of this departure is the absence of SPPs in diamond. Without the SPP scattering, 
the electrons in graphene lose a major energy relaxation mechanism and stay hotter than 
otherwise. Consequently, it is not unreasonable to expect a smaller drift velocity on a non- 
polar substrate (such as diamond) than on a polar counterpart.— Indeed, the results show 
that the average graphene electron energy on the diamond substrate is 0.44 eV at 30 kV/cm, 
while it is only 0.25 eV on Si02- 

To better examine the observed impact on electron drift velocities, the temperatures in 
the graphene film (T g ) and at the interface directly below (T s ) are provided in Fig. 2 as 
a function of applied bias. Both T g and T s increase with the field but their magnitudes 
vary widely depending on the substrate. Clearly, Si02 shows the extreme case of Joule 
heating with very high T g and T s that is consistent with the results of Fig. 1. Due to the 
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poor thermal conductivity Ksi0 2 > the excess heat emitted by graphene electrons cannot be 
efficiently channeled through the substrate. The resulting increase in temperature induces 
stronger electron-phonon scattering and subsequently degrades the drift velocity. As for 
/i-BN, the rise in T g and T s is far more modest despite the very low out-of-plane thermal 
conductivity Kh-BN(z), which is in fact about the same order of magnitude as K Si02 . The 
discrepancy comes from the in-plane thermal conductivity Kh-BN{x,y) that is about two 
orders of magnitude larger. Consequently, the transferred heat in /i-BN can easily spread in- 
plane unlike in SiC>2, utilizing a much wider thermal channel. Indeed, the 3D iso-temperature 
profile illustrates a laterally extended distribution near the interface — a sign of efficient heat 
removal. Two substrates with high thermal conductivities, diamond and SiC (unpassivated; 
not shown), show even smaller deviations from room temperature as expected. 

In the case of SiC-H, the moderate increase in T g (and the subsequent velocity decay) 
has a different origin. While thermal transport in the substrate is excellent (owing to a 
superior «sic), the heat conduction across the interface with a relatively large r gs provides 
the bottleneck. This point is clearly illustrated in Fig. 2(b) by the largest T g — T s among 
those plotted; the substrate surface temperature T s stays near 300 K in the graphene/SiC-H 
structure. Another interesting observation in Fig. 2(b) is that the absence of SPP interaction 
is visible from the result of diamond. The comparatively large T g — T s (over those of /i-BN 
and SiC^) indicates a greater disconnect between the graphene film and the substrate in 
terms of heat transfer. Since an identical r gs is assumed (except SiC-H), the absence of 
additional heat path via SPP emission appears to be the main reason for the increased 
T — T 

Finally, an attempt is made to compare the simulation results with the experimental 
studies available in the literature. Since the detailed heat dissipation (thus, the extent 
of Joule heating) depends on such factors as the exact geometry, etc., that not only vary 
from sample to sample but also are not fully characterized, it is rather difficult to have a 
meaningful one-to-one comparison at the quantitative level. For instance, when the lateral 
dimension of the graphene film is much larger than the thickness of substrate dielectric- 
heat transfer through the structure could simply be projected to a ID problem with the 
outcome potentially much different from that modeled in the current calculation. The 3D 
effect such as the lateral heat spread would be absent and /i-BN would behave more like SiC>2 
with pronounced Joule heating. At the same time, the additional structural details including 
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the location and the dimension of metal contacts could alter the heat dissipation pattern 
that are not included in the current model. Consequently, we treat the graphene sample 
size as an effective parameter that is adjusted to provide a good fit with the experimental 
data (specifically, those on Si0 2 ).-^ 

Figure 3 shows the comparison of v sat as a function of electron density n. The results 
clearly indicate the increased prominence of Joule heating at large values of n. This is obvious 
as more electrons mean lager heat generation per unit area, which leads to elevated lattice 
temperature and reduced v sat . Consistent with the results discussed above, Si02 suffers the 
biggest impact and then SiC-H is the next while the rest (BN, SiC, and diamond) remain 
largely unaffected in the considered range (n <4x 10 12 cm -2 ). One particularly interesting 
point to note is the slope of decay. The slope deduced from the simulation appears to 
become steeper in a close correlation with the rise in T g and reaches the dependence 
for Si02 that matches well with the experimental data from Refs. y and |4]. This result, 
however, cannot be explained when the Joule heating is excluded as evident from the figure 
(see the solid line). Consequently, it strongly indicates that the so-called l/y/n decay does 
not come from the SPP energy of the Si02 substrate as originally suggested.- Rather, it 
is a manifestation of Joule heating in its entirety including the influence of SPP scattering 
characteristics specific outcome and not a general rule. 

Considering the results thus far, Si02 may not be a desirable choice as a substrate ma- 
terial. As for diamond, the lack of SPP scattering appears to result in the drift velocities 
substantially smaller than other candidates (such as BN and SiC). However, it is also the 
most immune from the Joule heating degradation for its high thermal conductivity and may 
have an advantage at very high carrier densities (> 10 13 cm -2 ). While the performance 
of BN and SiC are generally comparable, SiC-H shows the sign of elevated temperatures 
in graphene. This could cause a significant concern in device breakdown characteristics in 
addition to the channel velocity reduction. Since the dangling bonds at the surface tend to 
be terminated in one form or another (by hydrogen or other specifies), the realistic structure 
involving the SiC substrate may be more like SiC-H than the ideal case without termination. 
Thus, it appears that /j-BN provides the best characteristics among the studied to interface 
with graphene in electronic applications. 
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SWAN-NRI programs. JMZ acknowledges support from NSF under the IR/D program. 
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Figure Captions 



Figure 1. (Color online) Drift velocities vs. electric field for graphene on different substrates. 
The graphene sample is assumed to be 1 /im x 0.5 /im with a carrier density of 1 x 10 12 cm -2 . 
The impurity density is 5 x 10 11 cm -2 . 

Figure 2. (Color online) (a) Graphene lattice temperature T g for different substrates, and 
(b) temperature difference T g — T s between the graphene lattice and the top surface of the 
substrate as a function of driving electric field. The conditions are the same as in Fig. 1. 

Figure 3. (Color online) Saturation velocity v sat vs. electron density in graphene on different 
substrates. In the calculations, it is assumed that the graphene film of 1 fim x 1 /im is under 
an electric field of 30 kV/cm. The impurity density is 5 x 10 11 cm -2 . The experimental data 
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are for the case of Si02 substrate. 
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FIG. 1: Li et al. 
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FIG. 2: Li et al. 
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FIG. 3: Li et al. 
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